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ABSTRACT 

OPS is an on-line system developed by M. Greenberger et. al. at 
Project MAC. The present work provides a powerful and simple way to 
perform numerical manipulations and calculations within OPS. The 
program package is called CALCULAID. 

A method of executing algebraic assignment statements, of which 
MAD and FORTRAN assignments are a subset, is provided. When this 
assignment-statement ability is coupled with other features of the 
OPS system, such as unconditional transfers, general conditionals, 
and array and functicm declarations, most of the ability of a compiler 
language is provided. Because Che programs written in OPS are executed 
interpretively, OPS-3 programs can be changed and re-run imediately, 
without being recoinplled. 

The other elements of CALCULAID are a program for creating multi- 
ple linear regression models, rank-ordering and counting data, and 
finding roots to polynomial equations in one unknown. 

The applications of CALCULAID to the analysis of a round-robin 
scheduling model and to a process-control problem are discussed, and 
conclusions regarding the suitability of running coii4>utatlonal programs 
in an interpretive mode are drawn. 
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Chapter I 
INTRODUCTION 

This work was performed within the framework provided by Project 
MAC, an experiment in multiple-access computer systems. Project MAC 
provides the means whereby several persons may use the conq>uter at 
what appears to them to be the same time, but which in reality is a 
single-server queue in which everyone lines up for service. A user 
places himself in the queue by Issuing a conraand to run either one of 
the system programs (such as logging in to the system or compiling a 
program) or one of his own programs, which may be anything he has 
previously written. 

Each user has space reserved for him in a disk file which he can 
access directly from his console. It is from this collection of programs 
in his own directory that he may select a program to be edited, compiled, 
or run. 

The great advantage of a system such as Project MAC's is that each 
user has a fairly direct line to the computer, but when he is not using 
the machine someone else may do so. If none of the "on-line" users has 
any programs to be run, the system supervisor recognizes this fact and 
runs programs that have been submitted in the traditional way for 
batch processing. Thus there is no penalty to the user if he wishes to 
take some time to think about what to do next. 



The Project MAC system, known as CTSS (Compatible Time Sharing 
System) [Ij, provides a wide spectrum of commands for program com- 
pilation and file manipulation. Since the Inception of the project 
in July, 1963 progress has been made In the area of Multiple Access 
Computers. With the introduction of a new file system in August, 
1965 users were given significantly better file-handling commands. 
The OPS system [2], operating within CTSS, gives the ability to 
structure coosnands and change programs as they are being run. 

The program "package" called CALCUIAID provides a convenient way 
to specify and execute requests for mathematical computation and man- 
ipulations. It runs entirely within the OPS-3 system. In addition 
to a general assignment -statement capability, such as the execution of 
the statement 

X-Y+A*(Z + 5) 

several more specific programs have been provided. These include a 
multiple linear or sinqple polynomial regression program, a solver of 
polynomial equations in one unknown, a count -by- interval program, and 
a program to rank-order a given sequence of nimibers. 

These programs can be used singly, to perform individual tasks, or 
they can be combined together in "programs", where each step of the 
program is one of these macro operations. The system can be used either 
as a very large and fast "desk calculator" or as a powerful programming 
device which allows simultaneous creation and execution of programs. 

A brief introduction to the OPS-3 system will be given, followed 
by a detailed description of the elements of CALCUIAID. 



The OPS System 

The reader interested In a detailed description of the OPS system 
is referred to the OPS-3 manual [2]. There he will find extensive 
explanations of all system programs and fuiketions. The brief Intro- 
duction to OPS given here serves as b«cfcground to this thesis. 

OPS extends the Project M/IC time-sharing system in the areas 
of command structuring and data handling. The OPS user finds it easy 
to Integrate his own coimaands into the system and write programs 
consisting of these comnands. Sach step of the program consists of one 
operator. This program of operators (called a compound operator or 
KOP) can then be used in the same way as an individual operator giving 
the user great flexibility in tailoring the command structure to his 
Individual needs. Steps that are performed repeatedly can be made 
into a KOP and called like a subroutine. 

A call to a KOP can be a small step, if the KOP has only a few 
simple operators, or a very large one, if the KOP contains calls to 
other KOP' 8 which call further KOP's, etc. In OPS, this hierarchical 
command structure can be built up on-line, with frequent checking of 
intermediate results. 

OPS provides a convenient facility for data storage and retrieval. 
It uses common storage as a data area that can be dynamically allocated 
by the user from the console, while the program is running. One can 
define cells, vectors, matrices, 3-dlmenslonal arrays and functions, 
in fixed, floating point, or alphabetic mode. Programs have been provided 
within the OPS system which allow symbolic reference to these data 
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Chapter II 
OPERATIONS WITHIN CALCUUID 

ASSIGNMENT STATEMENTS 

One conclusion Inpllclt within the Project MAC time-sharing 
system is that an individual command should initiate a significant 
amount of work, without putting an undue strain on the user. There 
has been considerable eEophasis on writing programs that are powerful 
and easy to use. 

The same policy has been followed in the construction of SET, the 
assignment- statement executor. SET executes the standard arithmetic 
operations +, -, *, and /. The special feature of SET is that operations 
can be specified for entire arrays. SET does a considerable amount of 
Interpretation and processing to determine indices and control loops » 
thereby relieving the programmer of his responsibilities in this area. 

Relieved of the tedious Job of index- and subscript-tending, the 
user has more freedom to think about his own problem instead of concerning 
himself with the "housekeeping" chores. The computer does drudgery with 
far less fatigue than do humans. 

An assignment statement is a line of the form 

SETA-X + Y + Z* SQRT.(C+D) 

The input line is analyzed, parsed into the proper groupings, and 
executed. The symbols on the right-hand side are looked up in the OPS 
symbol table, the Indicated operations are performed, and the result is 
stored in the location A (also looked up in the symbol table) . 

11 
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first row and second column of A. One very useful property of SET is 
that if any of the variables mentioned in an assignment statement are 
arrays » the statement is executed for all possible elements in the arrays. 
For Instance, if B is a vector of 10 elements, the statement 

SET B - 
sets all 10 cells of B to zero. If A is a vector of 10 elements, the 
statement 

SET A - B 
places the first element of B in the first cell of A, the second element of B 
in the second cell, etc. 

If, however, A is only 5 cells long, the process stops after only 
5 transfers have been made. In general, the range of operation is limited 
by the smallest array mentioned. Single cells and explicit nuaibera are 
expanded as far as is necessary, as in 

SET B - 2 

which sets all the cells of B to the value 2. 

SET B - DELTA 

fills B with whatever value DELTA has when the statement is executed. 

The same rules apply to 2- and 3-dimensional arrays. If A is a 
5 X 10 matrix and B is a 10 x 5 matrix, only the overlapping 5x5 part 
is affected by the statement 

SET A - B 
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As soon as the first row of B is exhausted SET begins to transfer the 
second row. After five rows have been thus transferred no more rows 
are available In A and the process Is terminated. 

Rows and coltinms of matrices may be treated as vectors. A(l) 
means the first row of A, while A(0,1) means the first column of A. 
The reason for the subscript is that SET recognizes a subscript of 
to mean a subscript over which an Iteration is to be perforiMd, An 
omitted subscript is also considered to be an iteration subscript. Thus 
A(I) and A(I,0) mean exactly the same in a SET statement. The 0^*^ cell 
of an array is referenced by setting a variable, say X, to zero and 
then referencing A(X). 

The only limitation placed on subscripts is that a subscripted 
variable may not be used within a subscript. A(B4<:/D) is allowable 
(if B, C, and D are cells), but A(I(1)) is not. The limitation Is 
necessary because of the way in which iterations are handled. Any 
other kind of computation, including calls to functions, may be Included 
In a subscript. 

It is often convenient to have a vector of the positive Integers, 
for use in computing weights or averages. An implicit vector of Integers 
called INTEGR has been Included within SET, and it may be used like 
any other vector In a SET statement. If A is a vector of 10 elements, 
the statement 

SET A » INTEGR 

will put the numbers 1,2,3,4, ... ,10 into A. Modes (integer or floating 
point) are again handled properly by SET, 
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Staple Arithmetic Stateneote 

The four arithmetic operators -f, -, *, and /, laeaning addition, 
subtraction, multiplication, and division reapectively, atid .P., 
aeanlng raise to the po«er are allowed. For exanple, the stateioent 

SET A - (ftK) * (D-E,P.2)/5 

would be executed as follows: If 8-1, C-4, I>»3, and E»6 

A « (144) * (3-<6)^/5 
- 5*(3-36)/5 
» 3-36 

" -33 

The nuafijer -33 would be placed in A, SET performs all necessary mode 
conversions. If some of the syafeols on the right hapf^en to be stored 
as integers, they are converted to floating-point nuaijers before any 
calcuUtioos arc performed. If A is declared to be integer mode In the 
symbol table, the result will be converted to Integer form before it is 
stored . 

The MAD {3] convention of bracketing all operators (with the 
exception of +, -, *, and /) with periods has been followed. Thus the 
power operator, as already mentioned, is ",?.". Bracketing makes opera - 
ora easily recognizable, and allows simple extension of the operator 
vocabulary without using special characters. 

Arrays 

Subscripts are handled just as In FORTRAN and MAD. a(1) indicates 
the first element of vector A. A(l,2) Indicates the element found in the 
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Function Call* 

Calls can ba tamim to all functions axcapt thoia that ratum • 
result through the calling aequanea. The reason a inunction cannot 
change its arguisents Is chat it is called with polntera to copies of 
all its arguments, not with pointers to the argiatmts theawelves. 
Therefore when a function attempts to chance soewthing in its calliiig 
sequence, it changes the copy» not the origl6*l cell with which it was 
called. 

Calls are made as in HAD, with a period between the naiie of the 
function and the parentheses enclosing the atganents. Argwiiiints may be 
computed, and the nooiber of arguamsts is li«ited to 10. X£ 700 with to 
use a function that has not yet been loaded, type 

LQADO ALPHA 

where ALFBA is the naaw of the file containiiME the daalred »«4Mroutine. Then 

ENTERS FUHCT BBXA I 

where BETA is the entry point to the subroutine (it «ay be ^e eaaie as ALPHA) 
and I is the nuaber of arguiaents the function expects. DataiUd descriptions 
of these loading and declaration statMients are ceotaiiMtd in the 0?8 aanual [2]. 
Kow SET has all the information it needs, and a statement of the form 

SET A - BETA.(C,IH€<3),X) 

will be properly executed. A call to a function may nc^ contain another 
function call. For exai^le, 

SET A " SQRT.(C+5IM.(B» 

is not allowed. 
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Array Operations 

SET will execute a statement for all elements contained in argu- 
ment arrays. It can also perform matrix multiplication, take first 
differences of vectors > or "coiiq)res8" one vector with another. This 
gives the user mathematical power beyond that offered by a compiler 
language. A single statement can do what would take 5 statements in 
HAD, and the chances for error are reduced, especially since the 
progranmer can be unconcerned with the dimensionality of arrays. If 
he specifies the multiplication of one matrix by another, SET will 
check to see that the nuoiber of columns in the first agrees with the 
number of rows in the second. If they do not agree, an appropriate 
message is printed. 

Literal Vectors 

If a vector Is to be used only once, say for a logical compression, 
it is bothersome to have to enter it In the symbol table and then do a 
"TYPE IN*' or a series of SET's. Literal vectors allow a vector to be 
typed directly Into a SET statement. To enter a set of values into 
VECO. type 

SET VECO - (2,4,A,G4-D,VEC1(5)) 

Any cooqiutations within this "literal vector" are allowable, except that 
none of the elements can be of dimension greater than 0. Thus 

SET VECO - (A+B*SIN.(5),B,C.L0G.(2),VEC1,M,VEC2) 

is legal, but 
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SET VECO = (A,B,C,VEC1+VEC2,D) 
Is not. 

These console vectors may be used nearly anywhere in a SET state- 
ment. The restrictions are noted below. 

Infix Operations 

Matrix multiplication is best explained by example. Assume 

A : cell 

VECO: vector of 5 elements 

VECl: vector of 10 elements 

VEC2: vector of 10 elements 

VEC3: vector of 20 elements 

MATl: matrix 5 x 10 

MAT2: matrix 10 x 10 

MAT3: matrix 10 x 10 

MAT4: matrix 10 x 20 

The statement 

SET A = VEC1.M.VEC2 

computes the inner product of VECl and VEC2 and puts the result in 
location A. 

SET A = VEC1.M.VEC3 

is an error since VECl and VEC3 are not of the same dimension. Constants 
may be used in matrix multiplications. 

SET A = VECl.M.l 
would put the sum of all elements of VECl into A. One possible application 
for this statement is the computation of sum-of -squares. One would use a 
sequence of statements like 
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SET VEC2 - VEC1*VEC1 
SET A - VEC2.M.1/10 

The first statement squares the elements of VECl, and the second sums 
them and divides to obtain the average squared deviation. These two 
statements could not be written 

SET A - (VEC1*VEC1).M.1/10 

because SET does not have the temporary storage that Is required. If 
the RMS were desired we could either add the statement 

SET A - SQRT.(A) 

or change the second statement to read 

SET A » SQRT.(VEC2.M.1/10) 
Matrices can be pre- or post -multiplied by vectors of the proper 
dimension. 

SET VECl - VEC2.M.MAT2 

performs the multiplication Indicated by 



VEC2 



MAT2 
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SET VECl - MAT2.M.VEC2 
perforins the multiplication 



MAT2 




The statement 

SET VECO » VEC2.M.MAT1 

Is allowable, since VEC2 and MATl are compatible, but the operation 
will be carried out for only the first 5 elements of VEC2 and the first 
5 columns of M^Tl. At that time there is no more space in VECO and the 
statement is terminated. 

The order of arguments in matrix-algebra is Important. 

VECO.M.MATl 

is allowable, since we are performing 

/ VECO 



MATl 
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MATl.M.VECO 
is not allouable, since it implies 




The 10 elements in each row of MATl are incompatible with the 5 
elements in VECO. 

Matrix multiplication by a constant is handled in a special way. 
If the left side of the SET statement is a vector » the constant is 
considered to be a vector also; if the left side is a single cell, the 
matrix is summed over all its elements and multiplied by the constant. 
For example 

SET VECl = MAT2.M.1 

performs 



HAT2 
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That is, VECl is set to the rofw sums of MAT2. 

SET VECl * 1.M.MAT2 

places the column sums of KA.T2 into VECl. The number 1 Is expanded as 
far as is necessary to satisfy MAT2. 
The statement 

SET A = 2.M.MAT1 

multiplies each element of MATl by 2 and cumulates the total into A. 
This seems the most legitimate interpretation of an operation which Is 
not defined in standard matrix algebra. 

NOTE: There is no possibility of confusing the period associated with 
an operator with a decimal point attached to a number. 

1.5.M.VEC1 

will be interpreted correctly, as will 

VECl. M.. 5 

Since the matrix multiplication operator automatically transposes a 
vector if it is required, there is no particular need to transpose 
vectors explicitly. .T. , which will transpose a matrix, has been in- 
cluded within SET. 

SET MAT2 = .T.MAT3+5 

will transpose MAT3, add 5 to each element, and put the result in MAT2. 
If .T. is given a non-matrix operand, SET will deliver an appropriate 
message. 



22 



An implicitly transposed matrix may not be used as an operand 
to ".M." 

SET VECl - (.T.MATD.M.VECO 

is not allowed. To perform the calculation, define a new symbol, 
say MATS, of dimension 10 x 5. Then type 

SET MATS - .T.MATl 
SET VECl - MATS;.M.VECO 

First-differences can be taken with the .D. operator. It has a 
single operand, which must be the name of an internally- stored vector 
(again, not a literal vector). It returns the series of first differ- 
ences between elements. Suppose VECl • .5,3.4,2.9,13.1,12.9,15.7,9.3. 
The statement 

SET VEC2 =- .D.VECl 

would place 2.9,- .5, 10. 2, -.2, 2. 8, -6. 4 in VEC2. Note that the vector 
generated by .D. has one less element than its operand. 

Logical Operations 

SET can perform the logical con^jression specified by Iverson [4] , 
which uses one array as a mask and another as the operand. The operator 
is called .C. and is used in expressions of the form 

VEC1.C.VEC2 
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If an element of VECl is positive, the corresponding eleoent of VEC2 is 
Included; if the element of VEGl is zero or negative, that elenent in 
VEC2 is not Included. For exaiiq;>le, suppose VECl were 0,1,1,0,0,0,1 »1,1 
and VEC2 were 1,2,3,4,5,6,7,8,9,10. 

SET VECO » VEC1.C.VEC2 

sets VECO to 2,3,7,8,9. If VECO is less than 5 elements long, the 
statement terminates before all 5 selections are made. As with matrix 
multiplication, .C. will not accept computed arguments. 

Even If the coiiq>resslng vector has only one bon-sero element, SET 
expects to find a vector on the left-hand side of the •" sign. At the 
time SET is checking the dimensionality of arguments it has no infor- 
mation about the contents of the compressing vector. 

CCHPUT 

SET can be used as a desk calculator. CCHPUT followed by any 
algebraic expression that can appear to the right of an - sign In a 
SET statement causes the result of the computation to be printed. 

COMPUT 2*(Y-Z.P.2) 

will print the result of the computation. The number will not be stored 
in the coiiq>uter. 

SETBCD 

BCD Information may be entered into single cells by SETBCD. 

SETBCD A - CAT 

places the word CAT into the cell A. No operator except the » may 
appear in a SETBCD statement. 
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REGRESSION MODELS 

The most popular kinds of models are linear models since they are 
most tractable analytically even though only an approximation to the 
real world. A typical modeling approach involves theorizing about the 
behavior of the system and deciding which variables are relevant, 
gathering experimental data concerning the behavior of the variable, 
and seeing how closely the data fit the model. 

At this point, the researcher can decide the system is no longer 
interesting; he can theorize a different kind of model; he can gather 
more data and try again; he can try various subtle manipulations with 
the model he has just created; or he can decide the model is good 
enough and go on to further work. 

The regression program in CALCUIAID is intended to assist in the 
last steps of the initial experimentation and In the subtle manipulations 
in the evaluation stage. The program is easy to use, and allows the 
user to make a change in the model, run it, evaluate it, make another 
change, etc. He will continue the cycle of modification and evaluation 
until he is satisfied with the results. 

Data Storag e 

FIT (the regression program) assumes that its data are stored in the 
following fashion: each variable, dependent or independent, must be stored 
in its own vector named to have mnemonic value to the user. The first 
element in each vector should represent an observation taken at the same 
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time as the first elements of all the other variables. Thus each vector 
is a time series of a different variable. It is crucial that the vectors 
all be of the same dimension, and that there be no unused locations 
at the ends of the vectors. FIT has no way of knowing how many points 
to include in the analysis except by examing the dimensionality of the 
vectors. 

Model Specification 

Suppose you were interested in a system which had 10 interdependent 
variables, and you had gathered data appropriate to this analysis. For 
ease in writing, assume the names of the variables (and the names of the 
vectors in which they are stored) are A,B,C,...,J. Initially you might 
want to theorize that A is really the independent variable, and all the 
others are independent. However, B and C are the most important ones, 
as they seem to influence A most directly. Type in 

FIT A TO B C 

This will perform the following analysis: 

Assume A is related to B and C linearly. That is, A is governed by 
a relationship A=c B+C2C+c^. FIT analyzes the vectors A,B, and C, and 
chooses those values for c , which give the most consistently close 
approximation to A. The criterion for consistency is that the sum of the 
squared deviations of A (observed) from A (calculated) is a minimum. The 
results are typed out in equation form, along with the calculation of the 
error. After examining the results you may want to include more variables 
in the model to reduce the error, or try a different combination of 
variables. 
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If some of the data points are suspect, you may want either to omit 
them from the analysis, or reduce their role In the model. This is 
accomplished by defining a vector, say W, which is of the same dimension 
as the variable vectors. Each element of W should be set to indicate 
the relative role each observation should play in determining the 
coefficients of the model. This weighting vector is called into play 
by including the word WEIGHT in the list of variables and following it 
by the name of the weighting vector, as in 

FIT A TO B C WEIGHT W 

or 

FIT A TO B WEIGHT W C 

The placement of the word WEIGHT is not important, as long as it occurs 
after the word TO. The two analyses, with the weights and without them, 
can then be compared. 

A vector into which residuals will be placed can be specified in 
the same way as the weighting vector. The key word is RESIDS, followed 
by the name of a vector of at least the dimension of the variable vectors. 

FIT A TO B RESIDS C 

will store in C how much each observation of the dependent variable differs 
from the value predicted by the model. It would be easy to change the 
weight attached to a given point and re-examine the residuals to deter- 
mine how much the error has changed. 

When a FIT is performed the mean square error is printed. Also 
printed is an "unbiased" value of the error which takes account of the 
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fact that each additional variable removes one degree of freedom from 

the system. For example, if we had three data points and three variables, 

an equation of the form 

A=c^B+C2C+C- 

could be found which would fit the three points exactly, no matter 

what the observed values were. The unbiased estimate then is obtained 

by multiplying the original estimate by N/(N-V-1), where N is the number 

of data points and V Is the number of variables in the analysis. In the 

example above, N/(N-V-1) is undefined, since N-V-1=0. If we fit A to 

B alone, the original estimate of error would be multiplied by S/O-l-D^S. 

Simple polynomial regression is a special case of linear regression, 
where the independent variables are all powers of a single variable. The 
specification of polynomial regression is 

FIT A TO B DEGREE N 

The word DEGREE is recognized by FIT and can appear anywhere on the line 
after the word TO. N may be specified symbolically or literally, and 
vectors of weights and residuals may be specified. 

FIT A TO B WEIGHT W RESIDS R DEGREE 2 

is quite acceptable. A will be fitted to a second-degree polynomial 
in B, using weights W. The residuals will be placed in R. 
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ROOTS TO POLYNOMIAL EQUATIONS 

■nie operator HOOT finds a root of a polynomial equation in one 
unknown. The method is the Newton-Raphaon technique of confuting the 
tangent at a particular point and using that as a local approximation 
to the curve. 

The solution is iterative, and it is possible to have a non-ter- 
minating loop. To avoid this, the user is asked to specify stme raaximum 
n^asiber of iterations. He is also asked to specify the tolerance or 
precision of the solution. The precision is the difference between 
successive approxinations . 

The parameters of ROOT are as follows: 

1. The power, P, of the highest-order term in the equation. This 
need not be an integer. It is assumed that all succeeding 
powers differ from this oaa by 1. That is, the equation is of 
the form 

P P-1 P-2 
a^x + a^x + ax + . . . 

2. The name of the array in which the coefficients are stored. 
The coefficients of the highest powers are first, with the 
constant term last. For exaiqple, if the equation were 

5x^*^ + 3x^*^ + 7x^'^ - 8 - 

the coefficient array would contain 5,3,7,-8. The coefficient 
array must be properly dimensioned in the symbol table. In 
this case, the vector must contain 4 cells. 

3. The starting value. This is the value used by ROOT in 
starting the analysis. It could be selected from a graph 
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of the curve, or from aomm other prior fe«llj% about the 
nature of the roots. If there ere eeverei roots to the 
equation they can be detcmlned by trying different starting 
values. 

4. The tolerance. The analysis Is coi^lete lAen the absoli^e 
value of the difference between successive s^proxla4tl<ms is 
less than this nuad>er< 

5. The iMTfliff^ number of Iterations. If the tolerance criterion 
Is not satisfied before the maxlamtm niMber of iterations has 
been performed, ftOOT retisras a "fence" as the valtte; 

Solution of the equation in the exas^le mentloBed above is 
acco^llshed In the following nay. the statement 

SET OO^F • (5,3,7,-8) 

will store the coefficients In COCFF, aad executing 

ROOT 3.5 COEFF S9kL TOL MUCIT 

will find the root of the equation. The numbers stored la SVkL^ TOL, 
and HAXn will be used as starting value, tolerance, and maxlitta Iter- 
ations respectively. The last three variables msy be entered directly, 
as in 

ROOT 3.5 aWFF 7.5 .00001 10000 
or thay may be omitted entirely. In irfilch case SVAL^, TOL". 001, and 
HfcXIT«1000 will be used. 

PRINT ROOT 3.5 COBTV 
will print the value on the console. 
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RANKING AND COUNTING 

Two operators have been provided to generate Information about the 
relative magnitudes of experimental data. The first, called RANK, rank- 
orders elements of one vector, A, into a second vector, B. The largest 
element In A is the first element of B, etc. The mapping (which element 
of A was the largest, which second largest, etc.) may be specified as 
the third parameter of RANK. The three parameters are 

1. the original vector 

2. the resultant vector (* if not desired) 

3. the vector of the mapping required to obtain the 
second vector from the first. 

For example, if A contained 8,9,3,7,4,6,2, execution of the line 

RANK ABC 

would place 9,8,7,6,4,3,2 into B and 2,1,4,6,5,3,7 into C. 

The second operator, called COUNT, places into vector B information 
as to how many elements of vector A fall into intervals specified by a 
third vector C. The intervals may be specified at the console directly. 
For example, suppose A contains 3,100,18,25,16,73,22 and 

SORT ABC 

is executed. If C contains 0,25,40,80,200 the first four elements of B 
are set to 4,1,1,1. 3,18,16 and 22 all fall in the interval 0-25, 25 
falls in the interval 25-40, 75 falls in the interval 40-80, and 100 
falls in the interval 80-200. 
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If C contains 0,200,25, B is set to 4,1,0,1,1,0,0,0. If the third 
element of C is less than the second, C is interpreted as specifying 
a range and a spacing. In this case, intervals between and 200 with 
spacing of 25 have been specified. 

If any elements of A fall outside of the intervals specified by C, 
an appropriate message is typed on the console. 



Chapter III 
APPLICATIONS 

Analysis of Round-Robin Scheduling In Time-Sharing 

An analysis of round-robin scheduling of requests for the services 
of a multiple-access computer affords an example of the use of CALCULAID. 
The analysis [7] was successful and was carried out quickly and easily 
with CALCULAID. 

The Scheduling Algorithm 

The scheduling algorithm works as follows: a first-in, first-out 
queue of users Is maintained. When a user desires the services of the 
computer he is placed at the end of the queue. The computer services 
users by removing the user at the head of the queue for either a basic 
time unit, called a quantum or to completion of his request which ever 
is shorter. When a user's request la finished, he is removed from the 
queue. If his request is not finished the balance is placed at the end 
of the queue, and he must wait until every other person in the queue has 
been run before he gets another turn. 

The sizes of requests are assumed to be distributed exponentially, 
as is the distribution of think times (the time between the completion 
of a user's request and his initiation of a new one). These assuiiq)tion8 
are supported by Scherr [5]. 
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The Analysla 

Four KOP's, or coa^otmd operators, «r« used la the calculation. 
The fir at, SETUP, aake for infornation (ii^ut para»tttera) and doea 
preliminary calculations. It calls PRj(»8, which calculates a vector 
P where P(I) is the steady-state probability of having I users in the 
queue for co^>letlon of pending requests, ft la th« probability of the 
queue being eapty. 

The calculation of the P(I> is dona in the loop beginning on line 
30 and terminating on line 60. PZ is calculated oA llxie 80, and Z is 
set equal to the sum of PZ and P^ by the assignMSnt 

Z-PZtf.M.l 

Z is then used to nomalise P and PZ so t^^ the p^tobabllltles sum to 
one. The most likely nun^er of persima In the qtt«M c«a be determined 
by ranking P (line 110). The first eleai^t of HfcP la tbe index in P 
of the largest element of P. If PZ la Urger ^UOi VQttrCl)) then it 
is most likely that the queue will be empty. 

CYCLE calculates YZ, the steady-state wait of the uaer for his 
first quantum, and YL, the ataady-aCata wait fee «U au«c«adiag quanta. 
The waits are obtained by calculating the aqH3>^ of pe^te in the queue 
and multiplying by the average request slse. COSTS computes the averaga 
cost accrued in the system for each request. Hote the call to the 
e^onentlatlon routine on line 10. After the coats hav« been calculated, 
control returns to OPS. Another run may be made by calling SETUP again. 
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KOP SETUP 

10 SET I - , II - 

20 TEXT HOW MANY Q 

30 TYPE IK FLOAT NQ 

40 TEXT ENTER THE Q 

50 TYPE IN FLOAT Q 1 Q NQ 

60 TEXT HOW MANY GAMMA 

70 TYPE IN FLOAT NG 

80 TEXT ENTER THE GAMMA 

90 TYPE IN FLOAT GAMMA 1 GAMIA NG 

120 SET I » I + 1 

130 IF I ,G. NQ 

140 RETRNK 

150 CALLK PROBS 

160 CALLK CYCLE 

170 TEXT FOR A Q OF 

180 TYPE OUT FLOAT Q I 

230 SET II » II + 1 

240 IF II .G. NG 

250 GOTO 300 

260 CALLK COSTS 

280 GOTO 230 

300 TEXT SINGLE LEVEL COSTS ARE 

310 TYPE OUT FLOAT CE 1 CE NG 

320 SET II = 

330 GOTO 120 
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KOP PROBS 

20 SET J = 1 , P ( N ) - 1 

21 SET X » EXP . ( - SIGMA * Q ( I ) ) 
23 SET X = ( 1 - X ) / SIGMA 

25 SET SIGMAP »(1-X)/(S+V) 

30 SETP(N-J)=P(N+1-J)* SIGHAP / ( ALPH * J ) 

40 SET J = J + 1 

50 IF J .L. N 

60 GOTO 30 

80 SET PZ - P ( I ) * SIGMAP / ( ALPH * N ) , Z » PZ + P .M. 1 

90 SET P » P / Z , PZ - PZ / Z 

100 SET RT - N / ( SIGMAP * ( 1 - PZ ) ) - 1 / ALPH 

110 RANK P * MAP 

120 TEXT MOST LIKELY NUMBER IN QUEUE IS 

130 IF P MAP 1 .G. PZ 

140 GOTO 170 

150 SET LIKELY - 

160 GOTO 180 

170 SET INDEX - MAP(l), LIKELY - P(INDEX) 

180 TYPE OUT FIXED LIKELY 

190 RETRNK 
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KOP CYCLE 

10 SET Z - 1 - P ( 1 ) / ( 1 - PZ ) 

20 SET YZ = ( S - Q ( I ) * X ) / SIGMA +V*(S+V/2) 

30 SET YZ=Z*YZ/(S+V) 

40 SET LB = N / ( 1 - PZ ) - SIGMAP / ALPH - 1 

50 SET Yl = YZ + ( LB - Z ) * ( S + V ) 

70 SET WP = YZ / ( 1 - ( S + V ) * SIGMAP * ( 1 - PZ ) ) 

80 RETRNK 



KOP COSTS 

10 SET Z = EXP . ( ( SIGMA + GAMMA ( II ) ) * Q ( I ) ) 

20 SET CE ( II ) = Yl + LB * ( S + V ) / ( Z - 1 ) + 
V / ( 1 - 1 / Z ) 

30 RETRNK 
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Model of a Cheese-Production Process 

Food production is frequently not well understood, because there 
are often too many variables in a real-world process to allow an 
analytic solution. An approximation in the form of a model may still 
be useful, and CALCUIAID can help in the formulation of the model. 

The Process 

Cheese-making is an art, say those who make cheese, because there 
are so many intangibles. The milk used may vary in age or in butter- 
fat and solids content. There are variations in process temperature and 
the amount of yeast used may vary. Air temperature may affect the 
process. The dependent variable is the yield; that is, what percentage 
of the milk turns into cheese. 

The table below contains observations which might have been taken 
over a period of several months. We will synthesize from the data a 
model that will help to predict future yields. The model could then 
be used in a linear program to optimize yields. 

The variables are stored in vectors of 20 elements each. Age of 
the milk is in the vector AGE, butterfat content in FAT, solids content 
in SOLIDS, process temperature in TEMP, amount of yeast in YEAST, air 
temperature in AIR, and percentage yield in YIELD. It should be mentioned 
that these data are purely hypothetical and almost certainly do not 
reflect reality. 
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Age of 
Milk 


Percent 
Butter fat 


Percent 
Solids 


Process 
Temperature 


Asiount of 
Yeast 


Air 
Temperature 


Yield 


10 


4.3 


2.3 


75.0 


8.7 


72 


15.42 


6 


6.2 


3.3 


73.4 


8.6 


73 


13.74 


4 


7.1 


2.5 


74.0 


7.8 


85 


11.98 


8 


5.4 


3.4 


73.7 


7.9 


84 


14.43 


3 


5.5 


3.1 


76.3 


8.5 


90 


9.11 


7 


3.9 


2.7 


76.2 


8.5 


74 


11.54 


14 


4.4 


2.7 


74.3 


8.3 


77 


19.49 


11 


5.1 


2.4 


75.8 


8.4 


70 


16.80 


6 


4.7 


2.8 


75.6 


7.8 


65 


11.10 


7 


4.8 


3.2 


76.1 


8.1 


68 


12.23 


2 


5.9 


3.0 


74.7 


8.6 


74 


9.05 


13 


4.2 


2.4 


73.9 


8.4 


75 


18.47 


9 


4.8 


2.1 


75.2 


8.2 


80 


14.56 


8 


5.3 


2.4 


75.8 


7.7 


83 


13.58 


12 


5.5 


2.3 


73.3 


7.9 


81 


18.65 


13 


3.8 


2.8 


76.9 


8.0 


88 


16.93 


5 


4.0 


2.5 


74.0 


7.7 


77 


9.82 


9 


4.9 


3.6 


73.8 


8.4 


77 


14.90 


10 


4.7 


3.5 


74.6 


8.3 


70 


15.70 


7 


5.7 


2.6 


74.7 


8.5 


72 


13.79 
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In the illustration, user Input vill be represented as lower 
case, computer output as upper case, and author coomient in parenthe- 
sis. It shows both the versatility of OPS and the particular abilities 
of FIT. 

fit yield to age 

YIELD - 

.8744 AGE + 
6.8947 
RMS DEVIATION - .89716 
CORRBCTED DEVIATION m .92047 

OK 

(To see If age is an important variable, try fitting yield 
to butterfat content.) 

fit yield to fat 

YIELD - 

-.9720 FAT + 

18.9341 
RMS DEVUTION - 2.91177 
CORRECTED DEVIATION « 2.8741 

OK 

(The RMS deviation is considerably higher, indicating that 
age is more Important than fat content. To aee if fat has 
much of an effect, fit yield to both age and butterfat.) 
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fit yield to age fat 

YIELD - 

1.0305 AGE + 

1.2055 FAT + 

-.4255 
RMS DEVUTION - .31955 
CORRECTED DEVIATION - .33684 

OK 

(The corrected RMS deviation was reduced from .920 to .337, 
a factor of nearly 3. These two variables have a consider- 
able effect on yield. See if the solids contents has much 
effect.) 

fit yield to age fat solids 

YIELD - 

1.0339 AGE + 
1.2054 FAT + 
.1061 SOLIDS + 
-.7475 
RMS DEVIATION <= .31645 
CORRECTED DEVIATION - .34323 

OK 

(The corrected deviation has increased, indicating that 
yield is not a function of the solids content. Try 
another factor, say the amount of yeast present.) 
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fit yield to age fat yeast 

YIELD - 

1.0312 AGE + 
1.2117 FAT + 
.5541 YEAST + 
-5.0141 
RMS DEVIATION - .26652 
CORRECTED DEVIATION - .28909 

OK 

(The deviation has decreased somewhat, though not 
strikingly. Include air temperature instead o£ 
yeast.) 

fit yield to age fat air 

YIELD - 

1.0338 AGE + 

1.2364 FAT + 

-.0157 AIR + 
.5979 
RMS DEVIATION - .30300 
CORRECTED DEVIATION - .32865 

OK 

(Again the deviation has increased, indicating air 
temperature Is not crucial.) 

This analysis Indicates that the yield Is related most strongly to 
the age of the milk, Its butterfat content, and the amount of yeast 
used. The analysis could be carried further, trying tranfformatlons 
of variables, more variables, or weighting son» observations more than 
others . 



APPENDIX 1 
Internal Charactertatlca of SET 

The group of subroutines known as SET executes assignment state- 
ments In the following steps: 

I. Parsing of the input line into meaningful subparts 

II. Creation of a modified "Polish list" (changing of the string 
from infix to suffix notation) 

III. Evaluation of the Polish list 

The names of the subroutines which perform these functions are BRKUP, 
HWaST, and EVAL. The routine which calls the subroutines in the 
proper order is called SET. tet us consider the routines in order. 

BRKUP 

BRKUP (short for BREAK UP) takes the input statement and delivers 
a line which is broken up into meaningful strings. The break characters 
+, -» *, /, (, ), ., », and , icoBom) are recognized. In addition, strings 
of less than 5 characters which are bracketed by "." are recognized as 
individual entities. A space is recognized as a break character but is 
not stored explicitly. 

The string 

X»A+SQRT.(C+D) 

would be broken up into 

X - A + SQRT . ( C + D ) 

42 
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The string 

VAR1'L0G.(B.P.3) 
would be broken into 

VARl = LOG . ( B .P. / 2 ) 
There is no possibility of confusing the . associated with an operator 
with the decimal point associated with a number. The string 

.P. .3 
is broken up into 

.P. .3 
and 

.P. 5 
is recognized as 

.P. 5 

MAKLST 

This subroutine accepts the broken-up string supplied by BRKUP and 
creates a modified Polish list. In the process it checks the symbol 
table for the dimensionality of the symbols appearing in the statement 
and assures that the statement is dimenslonally correct. 
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Polish Lists 

A Polish list is a prefix or suffix specification of a computation, 
and since SET uses suffix notation, that is what we will consider here. 
The list is constructed by assigning to each possible operator a 
precedence, which is an indication of its binding strength. For example, 
* (multiplication) has a higher precedence than +, and in the expression 
A+B*C the * will be performed before the +. 

The rules by which a Polish list is constructed from an inflx- 
notated expression are as follows: 

1. If the next element in the expression is a symbol, put it 
into the Polish list 

2. If the next element in the expression is an operator 

A. If the precedence of the operator is greater than or 

equal to the precedence of the last operator in the operator 
list, put the current operator into the operator list 

B. If the precedence of the operator is less than the 
precedence of the last operator, remove the last opera- 
tor from the stack and put it into the Polish list. Check 
the precedence again, and continue to remove operators 
until either condition A is satisfied or there are no 
more operators left. 

3. When the end of the expression is reached, put all operators 
into the Polish list. 

Consider the expression 

A+B*C-D 
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The first element encountered is A, a syinbol, so It Is put into the 
Polish list. +, with precedence 1, is put into the operator list. B 
is placed into the Polish list, and the status is 

Polish list A B 
Operator list + 

The operator * has precedence 2, which is greater than 1, so it is 
put into the operator list and C is put into the Polish list, giving 

Polish list ABC 
Operator list + * 

The - Is reached, and we go to step 2 B. * is removed from the operator 
list and appended to the Polish list, and - is added to the end of the 
operator list. We have 

Polish list ABC* 
Operator list + - 

D is added to the Polish list, and all operators are put into the Polish 
list (step 3). The final status is 

Polish list A B C * D - + 

and the operator list is empty. The evaluation of this list will be 
discussed in the next section. 

Modifications 

Some of the features of SET and of the OPS system require that more 
information be kept in the Polish list, information concerning sub- 
scripts of arrays, literal vectors, and function calls. 



46 

Any symbol that is dimensioned in the Syad>ol Table receives special 
treatment 4 It is followed in the Pollah list by the word SUB, which 
indicates that everything between it and the word SUB. , which must appear 
later, is part of a subscript. Each subscript is followed by the word 
FIXED) meaning the subscript is not to be iterated, or the word VAM, 
indicating it is to be iterated. SUB. is followed by a word which tells 
how many iteration subscripts the symbol has. 

If A Is a vector, A(l) appears in the Polish list as 

A SUB 1 FIXED SUB. 1 

If A is a matrix, A(l) appears as 

A SUB 1 FIXED 1 VARY SUB. 1 

and the second subscript of A is marked for iteration (omitted or zero 
subscripts are assumed to be Iteration •ubseripti). 

A(B4€,5) appears as 

A SUB B C + FIXED 5 FIXED SUB. 1 

and this will be correctly evaluated by fiVAL. 

Literal vectors are handled a little differently, the word CNSL. 
means that everything between itself and a matching CSX... constitutes a 
literal vector. The CSL.. is followed by the nuaiber of elements in the 
vector. The literal vector 

(1.3,7,F) 

is 

CNSL. 1 3 7 P I VA&Y CSt.. 1 
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the 1 VARY pUt U ufiftd by SVAI in its IttrMioA ifulgjclnt. *ai tha 1 
following CSL^. iadlc«t«a th« literal vtetor tui« 1 ItitAtloti ttibtctlpt 
(as do all lltatal vectors). The vector 

(B»Of3*P,A(5)) 

would api^aar ai 

CNSI. BC3f*4-AStmS rtXBO SUB. 1 CftL. . 1 

Function arnunanta are braekatad hy WW, at eha baiitming and a . 
at eha and. 

S<3Ri:.(A) 

baeoMMta 

sqitt fm. A . 
In eha Poiiah Uat. 

tha ttodifiad Poliah Hat ia fchan dttUvairad to fch« «¥*l4tation 
prograii^ 



Tha proj^aai WrkL avalnataa eha tuditiad Hli%h Hat ^ra^a^ad by 
KfclCXJf . A daacription o£ tha avaluaei<Mi »£ a atamdard Hlit^ liat 
will ba foUowad by an axplanation of tha pa«ull*'itiaa of KVAI*. 

tha an^raaaion 

A + B * C - D 
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is converted into the Polish list 

A B C * D - + 

EVAL simply examines each element of the list in turn. If the 
element is a symbol or a number EVAL stores its value in an evaluation 
list. If the element is an operator, EVAL executes the operator. 

In the present example, EVAL would encounter A, B, and C and their 
values would be stored In the evaluation list. When EVAL finds *, the 
last two elements In the list are removed and multiplied together, and 
the product Is put back on the list. It now contains A and B*C. D 
Is put on the list, and when - is found D Is subtracted from B*C. 
The evaluation list contains A and B*C-D. + Is picked up and executed 
by EVAL, giving A+B*C-D. 

Storing of Results 

EVAL must recognize that the execution of - Is not like the execu- 
tion of other operators. The first symbol in the Polish list is 
looked up but Instead of putting its value In the evaluation list, its 
address in the computer is inserted. Then when the - is executed the 
result Is stored In that address. 

Indexing and Iteration 

An internal function within EVAL, called INIT, uses the VARY words 
in the Polish list as indicators of iteration subacripts. It works like 
a speedometer, iterating on the last subscript until the symbol table says 
some array dimension is exceeded. It then increments the next higher 
subscript by 1 and resets the lowest level Indices to 1. 
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Every subscript of a particular level is incremented at the same 
time, and if any variable size is exceeded, execution of that iteration 
is terminated. The stopping criterion for the entire statement is the 
inability to increment any level of subscript and still stay within the 
size limits of all variables. 

Special Operations 

Certain operations in SET, namely .M. and .T.,must have more infor- 
mation than any of the others, which need only two floating-point numbers 
for their proper execution. .T. must have available the name and the 
indices of the matrix it is transposing. This information is determined 
by examining the Polish list, and the symbol is re-evaluated with the 
indices reversed. The new value then replaces the old one in the 
evaluation list. 

Execution of .M. is more complicated than that of .T. .M. must 
make decisions concerning the legitimacy of arguments and the dimension 
of results, and must handle iterations. This task is complicated by the 
fact that some of the iterations are handled by INIT and some are handled 
by special coding within the .M. section of the program, .M. examines 
the Polish list to determine the names and present indices of its 
operands and then performs the necessary multiplying and summing. This 
may involve multiplying one vector by another, a vector by a row or 
column of a matrix, a constant by a vector, etc. 
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Suggested Improvements to SET 

SET uses no significant temporary storage and as a result all opera- 
tions having to do with arrays are extremely slow. If A, B, and C 
are arrays and 

SET A = B + C 

is executed, the entire statement must be executed once for every element 
in B and C. The operations must be analyzed for every iteration. 

It would be more efficient to analyze the addition operation once, 
and perform the computations on the arrays at that time. CTSS routines 
can provide the temporary storage space required. 

The disadvantage of using temporary storage is that execution of 
the individual operators becomes more conqtlicated; execution of an 
arithmetic operation would involve more than processing two number». 
The speed tradeoff between faster handling of arrays and slower 
execution of operations would have to be determined by writing and 
testing the new program. 

Using temporary storage will remove the present restriction that 
some operations (.D., .C. , and .M.) will not accept coi)q>uted arguments. 
These operations will no longer go to the Polish list to pick up 
arrays, but will work, like other operations, on items in the evalua- 
tion list. The difference is that in using temporary storage, the items 
may be entire arrays instead of being restricted to single cells or 
nunibers. 
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APPENDIX 2 
Suitability of SET for Extensive Computatloo 

SET Is an Interpretive program, and each time a statement Is 
executed it must be cooQ>letely analysed. In addition* the analysis 
concerned with array operations Is quite wasteful If arrays are not 
being referenced. For exaiiq>le, execution of 

SET A-0 

Involves the execution of several hundred machine instructions. A 
compiled program would require two instructions. 

It is difficult to say exactly what the execution- time ratio 
between a coaplled program and a progran of SET* a would be, but for a 
small sample of test programs the ratio was about 1:100. A coiqpiled 
HAD program that took one or two seconds would take two or three minutes 
with set's. In the CTSS time-sharing environment this is not at all 
desirable, because three minutes of processor time would involve about 
half an hour of waiting at the console. A wait of that duration would 
be defeating to the CTSS aim of high-frequency user interaction, as 
well as being very wasteful of machine time. 

SET continues to be quite useful as a fast calculator, performing 
extended computations on-line, and as a debugging tool. However, after 
it is determined that a program (KOF) is running correctly. It is 
advisable to change it to a compiled program. This can be done in two 
ways. 
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The most obvious technique is aianply to manually translate the 
KOP into any standard language for which a cotapller la available. One 
could write a program which would look like the KOP version, using the 
same symbol names, etc. All SET statements which referenced arrays 
would have to be written out explicitly as loops, with the acconqpanying 
high probability of error in indexing. It may even be hard to find out 
Just what the dimensions of the symbols are, as SET looks them up at 
execution time. 

An easier conversion technique is MAOKOP [6]. This program 
transforms a KOP Into a MAD program which can then be compiled and 
loaded just like any other operator. The advantages to MADKOP are 
several. First, the entire translation, compilation and loading takes 
but a few seconds of processor time, instead of perhaps an hour of 
human time. Second, MADKOP can examine the syo^ol table as it is 
writing the MAD program, and will program SET'a iterations automatically. 
This relieves the user from having to worry at all about indexing. 

MADKOP provides a gain of about 25:1 in exacution time over straight 
interpretation. The full gain of 100:1 is not realized because MADKOP 
does not write in general the most efficient possible programs. It is 
fast, however, and relieves the user of most of his indexing responsibilities. 
A KOP should be compiled early in the debugging stage. If a KOP is to 
be executed often enough that it will use 10 seconds of processor time, 
it is almost certain to be advantageous to coaq>ile It. The con^llation 
and loading would take about 5 seconds, and total execution time would 
be less than half a second. 



53 



References 

1. The Compatible Time-Sharing System, A Programaer ' a Guide . Second 

Edition, The MIT Computation Center, P. A. Criaman, Editor, 
The MIT Press, Cambridge, Massachusetts, 1965. 

2. Greenberger, M. , et. al. The OPS-3 System for On-Linc Computation 

and Simulation . The MIT Press, Cambridge, Massachusetts, 1965. 

3. Arden, Galler, and Graham. Michigan Algorithm Decoder (MAD) . 

University of Michigan Press, Ann Arbor, Michigan, 1963. 

4. Iverson, K. A Programming language . John Wiley and Sons, Inc., 

New York, 1962. 

5. Scherr, A. L. An Analysis of Time-Shared Computer Systems . Project 

MAC, MAC-TR-18 (Thesis), June 1965. 

6. Morris, James. "Interpretive Systems in On-Line Programming". 

Unpublished Master's Thesis, Alfred P. Sloan School of 
Industrial Management, MIT, September 1965. 

7. Greenberger, M. The Priority Problem . A forthcoming Project MAC 

Technical Report. 



Associated References 

Hellerman, H. , "Experimental Personalized Array Translator System", 
Comm. of the ACM 7 (July 196A) , 

Hamblln, C. L., "Translation to and from Polish Notation", The 
Computer Journal (October 1962), pp. 210-213. 

Samelson, K. and Bauer, P. L., "Sequential Formula Translation", 
Conm. of the ACM (February 1960), pp. 76-83. 



This empty page was substituted for a 
blank page in the original document 



CS-TR Scanning Project 

Document Control Form Date : /cN 11 I HS 

Report # Lc5-Tr,-^o 

Each of the following should be identified by a checkmark: 
Originating Department: 

□ Artificial Intellegence Laboratory (Al) 
^ Laboratory for Computer Science (LCS) 

Document Type: 

^ Technical Report (TR) D Technical Memo (TM) 
D Other: 

Document Information Number of pages: ^±[£tdinsF^^^ 

' Not to include DOD forms, printer Wstrucdons. etc... original pages m^^ 

Originals are: Intended to be printed as : 

D Single-sided or D Single-sided or 

A Double-sided M Double-sided 

Print type: 

□ Typewriter □ OffsetPress □ LaserPrint 
"W InkJet Printer □ Unknown □ Other:__ 

Check each if included with document: 

^ DOD Fomi n Funding Agent Form ^ Cover Page 

D Spine D Printers Notes D Photo negatives 

D Other: 

Page Data: 



Blank Pages(byp^,.numb«):. 



Photographs/Tonal Material (by p.b« «"«*«):. 



Other (noli dMCfiptiorWkwga numtMO: 

Description : Page Number 






Scanning Agent Signoff: 

Date Received: /!/ // /^r Date Scanned: / / ^ Z?^ Date Returned: I l'> lU 



UNCLASSIFIED 



Security Classification 



DOCUMENT CONTROL DATA - R&D 

(Smcmlly cUaalUcmtlon ol Utim. body ol JiutrJCl and indmmini mtnotmllon mumt bt pnlntd «**n tltt avorrll raporl Im cUmmlHad) 



I. ORIGINATING ACTIVITY (Cotporml* muthot) 

Massachusetts Institute of Technology 
Project MAC 



Za. REPORT SECURITY CLASSIFICATION 

UNCLASSIFIED 



Zb. GROUP 



S. REPORT TITLE 

CALCULAID: An On-line System for Algebraic Computation and Analysis 



4. DESCRIPTIVE NOTES (Tn>m al rmpatt and Inetulw dalmt) 

Master's Thesis, Sloan School of Management 



S. AUTHOR<S) (Lamtntmm, Ur*t nMBt, InltlaS) 

Wantman, Mayer Elihu 



«. REPORT DATE 

September 1965 



7«. TOTAL NO. OF PAGES 

54 



76. NO. OF REFS 
10 



««. CONTRACT OR GRANT NO. 

Office of Naval Research, Nonr-4102(01) 

b. PROJECT NO. 

Nr-048-189 



••. ORIGINATOR'S REPORT NUMBERtSI 

MAC-TR-2G (THESIS) 



9b. OTHER REPORT NO(S) (Arry oOiat laMnbatm that mrnir bm 
aamltnad tbtm rmport) 



10. AVAILABILITY/LIMITATION NOTICES 

Qualified requesters may obtain copies of this report from DDC. 



11. SUPPLEMENTARY NOTES 

None 



12. SPONSORING MILITARY ACTIVITY 



Advanced Research Projects Agency 
3D-200 Pentagon 
Washington, D. C. 20301 



la. ABSTRACT 

OPS is an on-line system developed by M. Greenberger et. al. at Project MAC. 
The present work provides a powerful and simple way to perform nvunerical manipulations 
and calculations within OPS. The program package is called CALCULAID, and provides 
a method of executing algebraic assignment statements, of which MAD and FORTRAN 
assignments are a subset. When this assignment-statement ability is coupled with 
other features of the OPS system, most of the ability of a compiler language is 
provided. Because the programs written in OPS are executed interpretively, OPS -3 
programs can be changed and re-run immediately, without being recompiled. 

The applications of CALCULAID to the analysis of a round-robin scheduling 
model and to a process-control problem are discussed, and conclusions are drawn 
regarding the suitability of running computational programs in an interpretive 
mode. 



U. KEY WORDS 



Computer 

Machine-aided cognition 

Multiple-access computers 



On-line computer systems 
Real-time computer systems 
Time -sharing 



Time-shared computer systems 



DD /.??! 1473 (M.I.T.) 



UNCLASSIFIED 



Security Classification 



Scanning Agent Identification Target 



Scanning of this document was supported in part by 
the Corporation for National Research Initiatives, 
using funds from the Advanced Research Projects 
Agency of the United states Government under 
Grant: MDA972-92-J1029, 



The scanning agent for this project was the 
Document Services department of the M.I.T 
Libraries. Technical support for this project was 
also provided by the M.I.T. Laboratory for 
Computer Sciences. 



Sc&ntiedE 



Date: I f g^ mi 



M J«T* Librairles 
Document Services 



darptrgi.wpw Rev. 9/94 



